Pioneer Anomaly: Evaluating Newly Recovered Data 
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The Pioneer 10/11 spacecraft yielded the most precise navigation in deep space to date. However, 
their radio-metric tracking data received from the distances between 20-70 astronomical units from 
the Sun consistently indicated the presence of a small, anomalous, Doppler frequency drift. The 
drift is a blue frequency shift that can be interpreted as a sunward acceleration of ap = (8.74 ± 
1.33) X 10~^° m/s'^ for each particular spacecraft. This signal has become known as the Pioneer 
anomaly; the nature of this anomaly remains unexplained. 

New Pioneer 10 and 11 radio-metric Doppler data recently became available. The much extended 
set of Pioneer Doppler data is the primary source for new upcoming investigation of the anomaly. 
We also have almost entire records of flight telemetry files received from the the Pioneers. Together 
with original project documentation and newly developed software tools, this additional information 
is now used to reconstruct the engineering history of both spacecraft. To that extent, a thermal 
model of the Pioneer vehicles is being developed to study possible contribution of thermal recoil 
force acting on the two spacecraft. In addition, to improve the accuracy of orbital reconstruction, 
we developed a new approach that uses actual flight telemetry data during trajectory analysis of 
radio-metric Doppler flies. The ultimate goal of these efforts is to investigate possible contributions 
of the thermal recoil force to the detected anomalous acceleration. 

PACS numbers: 44.40.-|-a,45.20.df,95.10.Cc,95.10.Eg,95.55.-n,95.55. Jz,95.55.Pe,95.75.-z 



I. INTRODUCTION 

The first spacecraft to leave the inner solar system 
[J, S H, 3, Pioneers 10 and 11 were designed to con- 
duct an exploration of the interplanetary medium beyond 
the orbit of Mars and perform close-up observations of 
Jupiter during the 1972-73 Jovian opportunities. 

The spacecraft were launched in March 1972 (Pio- 
neer 10) and April 1973 (Pioneer 11) on top of identical 
three-stage Atlas-Centaur launch vehicles. After passing 
through the asteroid belt, Pioneer 10 reached Jupiter in 
December 1973. The trajectory of its sister craft. Pioneer 
11, in addition to visiting Jupiter in 1974^_also included 
an encounter with Saturn in 1979 (see [3, Q for more 
details). 

After the planetary encounters and successful comple- 
tion of their primary missions, both Pioneers continued 
to explore the outer solar system. Due to their excellent 
health and navigational capabilities, the Pioneers were 
used to search for trans-Neptunian objects and to estab- 
lish limits on the presence of low-frequency gravitational 
radiation 

Eventually, Pioneer 10 became the first man-made ob- 
ject to leave the solar system, with its official mission 
ending in March 1997. Since then, NASA's Deep Space 
Network (DSN) made occasional contact with the space- 
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craft. The last successful communication from Pioneer 

10 was received by the DSN on 27 April 2002. Pioneer 

11 sent its last coherent Doppler data in October 1990; 
the last scientific observations were returned by Pioneer 
11 in September 1995. 

The orbits of Pioneers 10 and 11 were reconstructed 
based primarily on radio-metric (Doppler) tracking data. 
The reconstruction between heliocentric distances of 20- 
70 AU yielded a persistent small discrepancy between 
observed and computed values [1, 0, [1] . After accounting 
for known systematic effects Q , the unmodeled change in 
the Doppler residual for Pioneer 10 and 11 is equivalent 
to an approximately sunward constant acceleration of 

ap = (8.74 ± 1.33) x 10"^° m/s. 

The magnitude of this effect remains approximately con- 
stant within the 3 dB gain bandwidth of the high-gain 
antenna (HGA) The nature of this anomalous ac- 

celeration remains unexplained; this signal has become 
known as the Pioneer anomaly. 

There were numerous attempts in recent years to pro- 
vide an explanation for the anomalous acceleration of 
Pioneers 10 and 11. These can be broadly categorized 
as either invoking conventional mechanisms or utilizing 
principles of "new physics" . 

Initial efforts to explain the Pioneer anomaly focused 
on the possibility of on-board systematic forces. While 
these cannot be conclusively excluded [1, 0] , the evidence 
to date does not support these mechanisms: the magni- 
tude of the anomaly exceeds the acceleration that these 
mechanisms would likely produce, and the temporal evo- 
lution of the anomaly differs from that which one would 
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expect, for instance, if the anomaly were due to thermal 
radiation of a decaying nuclear power source. 

Conventional mechanisms external to the spacecraft 
were also considered. First among these was the pos- 
sibility that the anomaly may be due to perturbations 
of the spacecrafts' orbits by as yet unknown objects in 
the Kuiper belt. Another possibility is that dust in the 
solar system may exert a drag force, or it may cause a fre- 
quency shift, proportional to distance, in the radio signal. 
These proposals could not produce a model that is consis- 
tent with the known properties of the Pioneer anomaly, 
and may also be in contradiction with the known prop- 
erties of planetary orbits. 

The value of the Pioneer anomaly happens to be ap- 
proximately cHq, where c is the speed of light and Hq 
is the Hubble constant at the present epoch. Attempts 
were made to exploit this numerical coincidence to pro- 
vide a cosmological explanation for the anomaly, but it 
has been demonstrated that not only this approach would 
produce an effect with the opposite sign [1, Q , it will also 
will have a much smaller magnitude. 

As the search for a conventional explanation for the 
anomaly appeared unsuccessful, this provided a motiva- 
tion to seek an explanation in "new physics". No such 
attempt to date produced a clearly viable mechanism for 
the anomaly [H. 

The inability to explain the anomalous behavior of 
the Pioneers with conventional physics has resulted in 
a growing discussion about the origin of the detected sig- 
nal. The limited size of the previously analyzed data 
set, also limits our current knowledge of the anomaly. 
We emphasize that in order to determine the origin of 
ap and especially before any serious discussion of new 
physics can take place, one must analyze the entire set of 
radio-metric Doppler data received from the Pioneers. 

In this paper we report on the progress of the recov- 
ery of the Pioneers' radio-metric Doppler data and flight 
telemetry and the status of the analysis of the Pioneer 
anomaly. The paper is organized as follows. In Section HIl 
we discuss the recovery of the extended Doppler data set 
and its current status. In Section Hill we present the flight 
telemetry and discuss its value for the upcoming inves- 
tigation. In Section HVl we discuss modeling the thermal 
recoil force and present preliminary results of this effort. 
In Section |V] we conclude with a summary and outline 
the next sets in the study of the Pioneer anomaly. 



II. RECOVERY OF THE EXTENDED 
DOPPLER DATA SET 

As of October 2007, an effort to recover all archived 
Pioneer 10 and 11 Doppler data, initiated at JPL in June 
2005, has been completed; there is now almost 30 years 
of Pioneer 10 and 20 years of Pioneer 11 Doppler data, 
most of which was never used in the investigation of the 
anomaly (Table |T|. 

The recovery process presented many unanticipated 



TABLE I: Pioneer 10 and 11 radio-metric Doppler data used 
in previous studies and presently available for new analysis. 



Data used in the previous analyses 


Spacecraft 
Pioneer 10 
Pioneer 11 


Time span 
03.01.87 - 22.07.98 
05.01.87 - 01.10.90 


Distances (AU) 
40.0 - 70.5 
22.4 - 31.7 


Currently available data 


Spacecraft 
Pioneer 10 
Pioneer 11 


Time span 
08.09.73 - 27.04.02 
10.04.73 - 11.10.94 


Distances (AU) 
4.56 - 80.2 
1.01 - 41.7 



challenges. Uncertainty in spin calibration and missing 
ramp information are two issues that created a major set- 
back in the data recovery project, as most of the newly 
recovered data files were not properly conditioned. In the 
following sections, we address these issues and briefly de- 
scribe the techniques used for orbital analysis. 



A. Spin calibration 

Radio communication with the Pioneer spacecraft 
was maintained using S-band microwave transmissions. 
The transmission frequency was ~2.1 GHz. The fre- 
quency of the return transmission from the spacecraft 
was ^2.3 GHz. When the spacecraft's radio communi- 
cation subsystem was in coherent mode, its receiver and 
transmitter were coupled, and the ratio of the frequency 
of the signal transmitted vs. the signal received by the 
spacecraft was exactly 240/221. It is this phase-coherent 
mode of operation that allowed precision Doppler mea- 
surements by ground stations, often with mHz accuracy. 

At this accuracy, one must take into account the com- 
bined effect of the spacecraft's spin and the polarization 
of the radio signal. 

The radio signal transmitted to, and returned by, the 
Pioneer spacecraft was circularly polarized. The space- 
craft was physically oriented such that its high-gain an- 
tenna always pointed towards the Earth; the spacecraft 
was spinning around its spin axis, which approximately 
coincided with the antenna axis. 

As a result, each revolution of the spacecraft added a 
full cycle to the signal in each direction. 

At a nominal spin rate of 4.5 revolutions per minute 
(rpm), this resulted a shift of 0.15 cycles per second in 
the return frequency. Knowing the actual spin rate of the 
spacecraft from telemetry, it is possible to calculate, and 
account for, the spacecrafts' spin more accurately. When 
a data file containing radio-metric information has been 
modified to account for this effect due to spin, it is said 
to be spin calibrated. 

When we assembled historical Pioneer Doppler data 
files for the first time, alarmingly we found that some 
of the files were spin-calibrated, whereas other files were 
not. It soon became evident that, in order to render 
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these files suitable for precision orbit calculations, their 
spin calibration had to be reconstructed from historical 
records (see discussion in Sec. IIIIPp . 



B. Ramp Data 

Several years after the launch of the twin Pioneer 
spacecraft, the Deep Space Network (DSN) began exper- 
imenting with "ramping" its transmission frequencies, in 
order to follow the variable frequency shift that is a result 
of the changing relative velocities of a distant spacecraft 
and a station on the Earth's surface, during the course 
of a data transmission session. 

Ramping is characterized by a starting frequency and 
a ramp rate, in Hz/s. Unfortunately, older data file for- 
mats used by the DSN had no provisions for storing this 
ramp data, therefore this information was stored sepa- 
rately. Further compounding the problem was the fact 
that ramp data records were often used to store the trans- 
mission frequency even when no actual ramping took 
place (i.e., the ramp rate was 0); the transmission fre- 
quency recorded in the Doppler record may or may not 
have coincided with the actual transmission frequency. 

When the recently recovered Doppler records were first 
assembled and analyzed, a strange "banding" appeared 
in the data, with data points offset by several hundred 
to several thousands Hz from each other. Missing ramp 
records were identified as a likely cause for this effect in 
many cases. In order to reconstruct a usable data set, it 
was necessary to search historical records that contained 
the necessary ramping information, or discard records for 
which no reliable ramp records could be found. 



C. Data Conditioning 

We realized that unusual behavior of the recovered 
raw Doppler data discussed in the previous section was 
due to the fact that these files were not submitted to 
a proper radio-metric data conditioning (RMDC) proce- 
dure [3| . This procedure is a standard part of preparing 
orbit tracking data for a trajectory analysis. As part of 
the standard process RMDC it includes both spin cali- 
bration and application of the correct ramp information 
to raw Doppler data received from DSN stations. 

It turned out that RMDC work was necessary for all 
Pioneer Archival Tracking Data Files (ATDFs; see dis- 
cussion in [5] ) that we received from National Space Sci- 
ence Data Center (NSSDC), as all of these files were raw 
unconditioned ATDFs. 

Below is a description of the RMDC process required 
for making Pioneer ATDFs usable for our investigation. 

The process is very detail-oriented, tedious, and labor- 
intensive. It consists of two phases - analysis and imple- 
mentation - that are discussed below. 




1970 1973 1976 1979 1982 1985 1988 1991 1994 1997 2000 2003 



FIG. 1: RTG 1 fin root temperatures (°F) for Pioneer 10. 



1. Analysis 

Analysis of tracking data files comprises the following 
critical steps: 

• Generate a summary of the contents of each ATDF; 

• Separate and print the tracking data, ramp data, 
and transmitter parameters of each file summary; 

• Compare ramp data to transmitter summaries to 
identify erroneous/extraneous transmitter on peri- 
ods, and uplink periods which are not covered by 
ramp records, then transfer this information to the 
data summary; 

• Check data summaries for erroneous reference fre- 
quencies (receiver and exciter) ; times when the ref- 
erence frequencies are being ramped, but there are 
no corresponding ramp records; erroneous status 
(validity, transmitter on/off, etc.) flags; erroneous 
identifier (exciter band, receiver band, Doppler 
mode, synthesizer, exciter type, receiver type, etc.) 
flags; erroneous Doppler bias values; 

Sometimes, the validation process can be an almost 
point-by-point exercise. 



2. Implementation 

After the analysis of tracking data files has been com- 
pleted, the necessary processing steps are implemented 
as follows: 

• Translate the analysis information into change 
commands (some of the necessary fixes could re- 
quire as many as 4 commands to implement) to be 
used with the appropriate program; 

• Execute the program to apply the correction com- 
mands (at least 250-(- commands for each ATDF); 
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• Generate a summary of the updated ATDF; 

• Check the new summaries to confirm that the cor- 
rections have been properly applied, and no new 
problems have been introduced (or "come to the 
surface") by making the changes. 

Once the steps above are implemented, the tracking 
data files are ready to be used for generating standard 
format Orbit Data Files (ODFs) that can be submitted 
for further analysis with the Orbit Determination Pro- 
gram (ODP) software. 

After the two-stage process above is complete, another 
person takes a "conditioned" ODF and runs it through 
ODP to verify the quality of the file. Experience has 
shown that this step sometimes identifies further prob- 
lems, resulting in some additional time and effort being 
needed before a fully-usable tracking data file is avail- 
able. If the quality is satisfactory then the file is added 
to the collection of good files, otherwise the file is re- 
turned for further RMDC. The process repeats until all 
files are found acceptable. 

The work on data conditioning of nearly 600-1- Pioneer 
ATDFs is progressing well and will be completed by the 
end of 2007. The result of this analysis is that at long 
last, a usable data set is now available for analysis, cover- 
ing in particular the early parts of the Pioneer 10 and 11 
missions, for which no Doppler data records were avail- 
able in the past. This is especially important because 
there has been speculation that there may have been an 
"onset" of the anomalous acceleration of Pioneer 10 and 
11 when they reached approximately the orbit of Saturn 
(see discussion in Q). Is this onset a real effect, or just 
an artifact of incomplete past analyses? If the effect is 
real, was it caused by, for instance, incorrect solar pres- 
sure calibration, or does it represent a hint about the real 
cause of the Pioneer anomaly? Our upcoming analysis is 
aimed to answer these and other questions. 
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FIG. 2: Changes in total RTG electrical output (in W) on 
board Pioneer, as computed using the mission's on-board 
telemetry. 



general relativistic corrections, non-gravitational forces, 
and equations describing the propagation of the radio 
signals through the interplanetary medium, the gravita- 
tional field of the Sun and planets, and the Earth's at- 
mosphere and ionosphere. If needed, other forces may 
be incorporated into the model and used for hypothesis 
testing (see details in [3]). 

The aim of the orbital analysis is to determine the 
initial position and velocity (state vector) of the space- 
craft, as well as the values of additional, model-dependent 
parameters that yield a best match with observational 
data. In our case, the goal is to reduce the difference be- 
tween calculated and measured Doppler values (i.e., the 
Doppler residual) and then adjust the initial state vector 
and other parameters until the residual is minimized (for 
details on models, methods, and formats consult |3|,l7|,[8|). 



D. Orbital analysis with Doppler data 

Doppler data are the recorded measurements of the 
shift in frequency of the signal received from the space- 
craft, as a result of the combined relative motion of the 
spacecraft and the Earth station(s). The position of 
the transmitting station changes as a function of time 
as a result of the combined effects of the Earth's rota- 
tion, precession, and orbital motion. The position of the 
spacecraft can be calculated if its orbital parameters are 
known. After these calculations, the Doppler shift in 
frequency can be estimated and compared with the ob- 
served frequency shift (for more details on the analysis 
of radio-metric Doppler data and formats used for these 
purposes consult [sj, 0] and references therein) . 

In principle, an orbit determination algorithm employs 
a generalized equation of motion (for extensive details, 
see 0), incorporating not just the Newtonian gravita- 
tional effects of massive solar system bodies, but also 



III. FLIGHT TELEMETRY 

The anomalous acceleration of the Pioneer spacecraft is 
very small. Only ~65 W of anisotropically emitted ther- 
mal power is sufficient to produce an acceleration of com- 
parable magnitude. As the Pioneer spacecraft utilized 
radioisotope thermoelectric generators (RTGs) that pro- 
duced several kW of waste heat, an anisotropy as small 
as ~3% in the pattern of thermal radiation can produce 
acceleration of the needed magnitude. For this reason, it 
is very important to characterize the thermal emissions 
of the Pioneer spacecraft as accurately as possible. 

A unique feature of the current effort is the use of 
telemetry files documenting the thermal and electrical 
state of the spacecraft. This information was not avail- 
able previously; however, by May 2006, the telemetry 
files for the entire durations of both missions were re- 
covered, pre-processed and are ready for the upcoming 
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FIG. 3: A drawing of the Pioneer spacecraft. 



study. Both of the newly assembled data sets are pivotal 
to establishing the origin of the detected signal. 

All transmissions of both Pioneer spacecraft, including 
all engineering telemetry, were archived 5] in the form 
of files containing Master Data Records (MDRs). Orig- 
inally, MDRs were scheduled for limited retention. For- 
tunately, the Pioneers' mission records avoided this fate: 
with the exception of a few gaps in the data |5|] the en- 
tire mission record has been saved, comprising ~60,000 
data points for Pioneer 10, and ~50,000 data points for 
Pioneer 11, a total of '^35 GB. These recently recov- 
ered telemetry readings are important in reconstructing 
a complete history of the thermal, electrical, and propul- 
sion systems for both spacecraft. This, it is hoped, may 
in turn lead to a better determination of the spacecrafts' 
acceleration due to on-board systematic effects. 

Telemetry formats can be broadly categorized as sci- 
ence formats versus engineering formats. Telemetry 
words included both analog and digital values. Digital 
values were used to represent sensor states, switch states, 
counters, timers, and logic states. Analog readings, from 
sensors measuring temperatures, voltages, currents and 
more, were encoded using 6-bit words. This necessar- 
ily limited the sensor resolution and introduced a sig- 
nificant amount of quantization noise. Furthermore, the 



analog-to-digital conversion was not necessarily linear; 
prior to launch, analog sensors were calibrated using a 
fifth-order polynomial. Calibration ranges were also es- 
tablished; outside these ranges, the calibration polyno- 
mials arc known to yield nonsensical results. 

With the help of the information contained in these 
words, it is possible to reconstruct the history of RTG 
temperatures and power, radio beam power, electrically 
generated heat inside the spacecraft, spacecraft temper- 
atures, and propulsion system. 



A. RTG temperatures and power 

The exterior temperatures of the RTGs were measured 
by one sensor on each of the four RTGs: the so-called "fin 
root temperature" sensor. Figure [1] depicts the evolution 
of the RTG 1 fin root temperature for Pioneer 10. 

A best fit analysis confirms that the RTG tempera- 
ture indeed evolves in a manner consistent with the ra- 
dioactive decay of the nuclear fuel on board. The results 
for all the other RTGs on both spacecraft are similar, 
confirming that the RTGs were performing thermally in 
accordance with design expectations. 

RTG electrical power can be estimated using two sen- 



sor readings per RTG, measuring RTG current and volt- 
age, from which power can be computed by direct calcu- 
lation. 

All this electrical power is eventually converted to 
waste heat by the spacecrafts' instruments, with the ex- 
ception of power radiated away by transmitters. 



B. Electrically generated heat 

Whatever remains of electrical energy (Fig. [2]) after 
accounting for the power of the transmitted radio beam 
is converted to heat on-board. Some of it is converted to 
heat outside the spacecraft body. 

The Pioneer electrical system is designed to maximize 
the lifetime of the RTG thermocouples by ensuring that 
the current draw from the RTGs is always optimal. This 
means that power supplied by the RTGs may be more 
than that required for spacecraft operations. Excess elec- 
trical energy is absorbed by a shunt circuit that includes 
an externally mounted radiator plate. Especially early in 
the mission, when plenty of RTG power was still avail- 
able, this radiator plate was the most significant com- 
ponent external to the spacecraft body that radiated 
heat. A specific telemetry word tells us the shunt circuit 
current, from which the amount of power dissipated by 
the external radiator can be computed using the known 
ohmic resistance of the radiator plate. 

Other externally mounted components that consume 
electrical power are the Plasma Analyzer, the Cosmic 
Ray Telescope, and the Asteroid/Meteoroid Detector. 
Though these instruments' exact power consumption is 
not telemetered, we know their average power consump- 
tion from design documentation, and the telemetry bits 
tell us when these instruments were powered. 

Two additional external loads are the battery heater 
and the propellant line heaters. The power state of these 
loads is not telemetered. According to mission logs, the 
battery heater was commanded off on both spacecraft on 
12 May 1993 [|. 

Yet a further external load is the set of cables connect- 
ing the RTGs to the inverters. The resistance of these 
cables is known. Using the RTG current readings it is 
possible to accurately determine the amount of power 
dissipated by these cables in the form of heat. 

After accounting for all these external loads, whatever 
remains of the available electrical power on board is con- 
verted to heat inside the spacecraft. So long as the body 
of the spacecraft is in equilibrium with its surroundings, 
heat dissipated through its walls has to be equal to the 
heat generated inside. 



C. Compartment temperatures and thermal 
radiation 

As evident from Fig. [31 the appearance of the Pioneer 
spacecraft is dominated by the 2.74 m diameter high gain 
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FIG. 4: Bottom view of the Pioneer 10/11 vehicle, showing 
the louver system. A set of 12 2-blade louver assemblies cover 
the main compartment in a circular pattern; an additional 
two 3-blade assemblies cover the compartment with science 
instruments. 



antenna (HGA) . The spacecraft body, located behind the 
HGA, consists of a larger, regular hexagonal compart- 
ment housing the propellant tank and spacecraft elec- 
tronics; an adjacent, smaller compartment housed science 
instruments. The spacecraft body is covered by mul- 
tilayer thermal insulating blankets, except for a louver 
system located on the side opposite the HGA, which was 
activated by bimetallic springs to expel excess heat from 
the spacecraft. 

Each spacecraft was powered by four radioisotope ther- 
moelectric generators (RTGs) mounted in pairs at the 
end of two booms, approximately three meters in length, 
extended from two sides of the spacecraft body at an an- 
gle of 120°. A third boom, approximately 6 m long, held 
a magnetometer. 

The total (design) mass of the spacecraft was ~250 kg 
at launch, of which 27 kg was propellant @ . 

For the purposes of attitude control, the spacecraft 
were designed to spin at the nominal rate of 4.8 rpm. Six 
small monopropellant (hydrazine) thrusters, mounted in 
three thruster cluster assemblies, were used for spin cor- 
rection, attitude control, and trajectory correction ma- 
neuvers (see Fig. ini). 

The passive thermal control system consisted of a se- 
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FIG. 5: Location of thermal sensors in the instrument com- 
partment of Pioneer 10/11 6]. Temperature sensors are 
mounted at locations 1 to 6. 

ries of spring-activated louvers (see Fig. [4|) . The springs 
were bimetallic, and thermally (radiatively) coupled to 
the electronics platform beneath the louvers. The louver 
blades were highly reflective in the infrared. The assem- 
bly was designed so that the louvers fully open when tem- 
peratures reach 30° C, and fully close when temperatures 
drop below 5°C. 

As the total exterior area of the spacecraft body and 
the effective emissivity of the thermal blankets covering 
the spacecraft are known [To| . a crude initial calculation 
can be made determining the amount of heat radiated 
by the spacecraft walls and through the louver system. 
These estimates turn out to be consistent with pre-launch 
thermal vacuum chamber tests [ill, of louver sys- 
tem, and with temperature readings obtained from on- 
board temperature sensors. 

There are 6 platform temperature sensors (Fig. [5]) in- 
side the spacecraft body: 4 are located inside the main 
compartment, 2 sensors are in the science instrument 
compartment. The main compartment has a total of 12 
2-blade louver blade assemblies; the science compartment 
has 2 3-blade assemblies. 

The thermal vacuum chamber tests provide values for 
emitted thermal power per louver assembly as a function 
of the temperature of the electronics platform behind the 
louver. This allows us to estimate the amount of thermal 
power leaving the spacecraft body throu gh the louvers, 
as a function of platform temperatures [l3|, providing 
means to estimate the amount of heat radiated by the 
louver system. 



D. Spin 

The spin rate of the Pioneer spacecraft was measured 
independently by three instruments: a star sensor and 
two Sun sensors. Each of these instruments was capable 



of providing a roll reference pulse, which in turn was used 
to synchronize other spacecraft activities. 

The spin rate was telemetered to the ground using 3 6- 
bit words. The resulting 18-bit value measured the time 
of one full revolution of the spacecraft in units of 1/8192 
seconds using an on-board crystal oscillator as a reference 
clock. 

Unfortunately, the star sensor on board Pioneer 10 
stopped functioning shortly after Jupiter encounter. The 
Sun sensors, in turn, could only operate if the Sun 
was bright enough and appeared sufficiently far from 
the spacecraft spin axis. Consequently, the Sun sensors 
stopped providing a useful roll reference pulse after the 
spacecraft reached a heliocentric distance of ~35 astro- 
nomical units. 

The Pioneer support team devised an alternate pro- 
cedure, utilizing the Imaging Photo-Polarimeter (IPP) 
instrument and its abilit y t o image faint objects in a 
precisely timed manner [14| . Later, when there was 
no longer sufficient power available on board to operate 
this instrument, the spin rate was estimated using other 
means, such as relying on the previously established rate 
of spin change, and on the spin introducing a slight mod- 
ulation of the spacecraft's radio signal when the spin axis 
did not exactly coincide with the spacecraft-Earth line. 

The star sensor and Sun sensors remained fully opera- 
tional on Pioneer 11 throughout its mission. 

As a result, we have high accuracy spin information 
for both spacecraft for much of their missions. The spin 
behavior of the two spacecraft was not identical (Fig. [6]). 
The spin rate of Pioneer 10 was decreasing steadily, and 
does not appear to be affected by maneuvers. The spin 
rate of Pioneer 11 increased rapidly early in the mission 
as a result of a spin thruster malfunction, but it continued 
to increase afterwards. A close examination of the spin 
behavior reveals, however, that the increases in spin rate 
are discontinuous and coincide with maneuvers; between 
maneuvers, the spin rate was decreasing (Fig.[7|). 



E. Maneuvers 

Each of the Pioneer spacecraft was equipped with six 
monopropellant hydrazine thrusters, designed to perform 
three types of maneuvers: spin-up/spin-down maneuvers, 
delta-y maneuvers, and precession maneuvers. The first 
of these was used early in their missions to adjust the 
spin rate after the deployment of the RTG and magne- 
tometer booms. The second type of maneuver was used 
to correct the spacecrafts' trajectories, to achieve the de- 
sired planetary encounters with Jupiter and (for Pioneer 
11) Saturn. Precession maneuvers, in turn, were used 
regularly throughout the mission; their goal was to ad- 
just the spacecrafts' spin axis, to ensure that the HGA 
continuously points in the direction of the Earth. 

In theory, precession maneuvers do not alter the ve- 
locity of the spacecraft, only change their orientation. 
In practice, due to small differences between thrusters. 
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FIG. 6: On-board spin rate measurements (in RPM) for Pioneer 10 (left) and Pioneer 11 (right). The sun sensor used on 
Pioneer 10 for spin determination was temporarily disabled between November 1983 and July 1985, and was turned off in May 
1986. Continuing spot measurements of the spin rate were made using the Imaging Photo-Polarimeter (IPP) until 1993. The 
anomalous increase in Pioneer ll's spin rate early in the mission was due to a failed spin thruster. 




Pioneer 1 1 telemetered spin rate (RPM) 



FIG. 7: Zoomed plots of the spin rate of Pioneer 11. On the left, the interval examined in ref. [3| is shown; maneuvers are 
clearly visible, resulting in discrete jumps in the spin rate. The figure on the right focuses on the first half of 1987; the decrease 
in the spin rate when the spacecraft was undisturbed is clearly evident. 



uncertainties in the duration of firing pulses, thruster 
misalignment, and outgassing after maneuvers, a small, 
essentially random change in the spacecrafts' motion is 
expected. These changes can be modeled, for instance, 
as an instantaneous change in the line-of-sight velocity 
of the spacecraft. However, before one can do such mod- 
eling, the time when the maneuvers occurred must be 
known. 

Flight telemetry contains information about the num- 
ber of pulses fired by each of the six thrusters on board. 
The number of firing pulses recorded, and the thrusters 
participating in the maneuver unambiguously identify 
the type of maneuver that was performed. 

Although thruster pulse counts are monitored contin- 
uously, due to low data rates, the corresponding data 



words are telemetered as infrequently as once every 51.2 
minutes at the lowest telemetry data rate. Due to gaps 
in ground station coverage or missing data records, there 
are instances when an increase in the thruster pulse count 
is bracketed by two telemetry readings that are many 
hours, perhaps even more than a day apart. However, 
due to the smallness of the velocity change as a result 
of a precession maneuver, even such an uncertainty in 
the timing of the maneuver does not adversely affect the 
accuracy of orbital modeling. 



FIG. 8: A geometric model (left) of the Pioneer spacecraft, used for finite element analysis, and a photograph (right) of 
Pioneer 10 prior to launch. The geometric model accurately incorporates details such as the Medium Gain Antenna (MGA), 
the Asteroid-Meteoroid Detector, and the star sensor shade. Note that in the geometric model, the RTGs are shown in the 
extended position; in the photograph, the RTGs are stowed. 



IV. MODELING THE THERMAL RECOIL 
FORCE 

Even with all the recently recovered telemetry infor- 
mation, modeling the thermal recoil force on the Pioneer 
spacecraft remains a difficult task. The main reason for 
this can be summed up as follows: What we seek is not 
just the amount of heat emitted by the spacecraft (which 
we know with good precision), but the small anisotropy 
in the pattern of its thermal emissions, an inherently 
second-order effect. 

The traditional way of estimating thermal recoil forces 
is to build a detailed thermal model using, for instance, 
finite element analysis tools. We are also exploring an al- 
ternative approach that combines the use of flight teleme- 
try with precision orbit calculations, and provides a more 
direct and, it is hoped, more accurate estimate of the 
anisotropy of thermal radiation p^ . 

A. Finite element thermal modeling 

The radiative exchange of thermal energy between two 
surfaces is described by the equation 

[ dMdA,. (1) 

JA2 JAi f 

where Q is the radiative power emitted by the surface 
Ax that is absorbed by the surface A2, I is the intensity 
of radiation at surface element dAi, r is the distance 
between surface elements dAi and dA2, and 9i and 62 are 
the angles between the line connecting surface elements 
dAi and dA2 and their respective normals. 



The intensity of radiation is related to the radiant 
emittance q by the formula / = ^/tt. The radiant emit- 
tance, in turn, can be calculated from the surface tem- 
perature T and emissivity coefficient e using the Stefan- 
Boltzmann law: 

q - atT\ (2) 

where a ~ 5.67 x 10"** Js^^m-^K-^ is the Stefan- 
Boltzmann constant. 

In principle, these equations completely determine the 
thermal emissions of an object. The anisotropy of ther- 
mal emissions can be determined, for instance, by sur- 
rounding the object with an infinite (or, at least, suffi- 
ciently large) black sphere and calculating the amount of 
radiation absorbed by this sphere using Eq. ([1]). 

In order to carry out the calculations, however, one 
must have accurate knowledge of 

• the emitting object's detailed geometry; 

• the object's temperature at all points across its sur- 
face; 

• the infrared emissivity coefficient e at all points 
across the surface, and its dependance on the tem- 
perature. 

This information is not always readily available. In the 
case of Pioneer 10 and 11, only a few detailed drawings 
have been found from which the spacecraft's geometry 
could be estimated. Surface properties are not known 
very accurately, and very little is known about the ef- 
fects of aging in the deep space environment. Although 
the spacecraft had several temperature sensors on board. 



10 



0.2C (C319PLT 5) 
vs 16.1C 



I vs 1 



5C (C318 PLT4) 
5.4C 



3.2C (C302 PLT2 
vs 13.4C 




^2.9C (C320 PLT 6) 
vs 15.4C 



0.2C (C304 PLT 3) 
vs 12.4C 



11. 9C (C301 PLT 1) 
vs 13.8C 



FIG. 9: A "work-in-progress" temperature map of the outer surface of the Pioneer 10 spacecraft body, comparing temperatures 
calculated via a numerical finite element method vs. temperatures measured by platform temperature (PLT) sensors and 
telemetered. While agreement between calculated and telemetered temperatures is expected to improve as the model is being 
developed, discrepancies between these values illustrate the difficulties of creating a reliable temperature map using numerical 
methods. 



these were designed to measure interior temperatures; 
the external temperature distribution can only be calcu- 
lated by constructing a complex model that takes into 
account the internal structure and heat conductive prop- 
erties of the spacecraft. 

Despite these difficulties, a highly detailed mathemat- 
ical model of the spacecraft has been constructed suc- 
cessfully [l3|. This model is used to calculate radiative 
heat exchange between spacecraft surfaces, absorbed so- 
lar loads, heat flows, and predicted temperatures. The 
finite element model incorporates ^3000 nodes and 2600 
plate elements, using 3.4 million radiation conductors 
and '^TOOO linear conductors. 

Results from this model can be validated against the 
readings from the spacecraft's temperature sensors, as 
seen in Figs. [H and [TOl 

Flight telemetry can provide another important means 
of model verification. In addition to temperature sensors, 
flight telemetry also offers readings from which the ther- 
mal power generated on board can be calculated. The 
principle of energy conservation dictates that in a space- 
craft that's in steady state, the amount of thermal power 
generated must equal the amount of heat emitted by the 
spacecraft. Because of this, temperature and power read- 



ings can be viewed as redundant parameters characteriz- 
ing the same physical processes. 



B. Estimating thermal recoil forces from orbital 
analysis 

Our analysis of the geometry and thermal properties 
of the Pioneer spacecraft yielded the following flndings: 

• The thermal output of the Pioneer spacecraft can 
be accurately modeled using as few as two heat 
sources, electrical and RTG heat; 

• As the spacecraft is spinning, lateral forces result 
in no significant long-term acceleration; the spin 
axis, in turn, points approximately in the direction 
of the Earth at all times; 

• The acceleration of the spacecraft due to thermal 
radiation is a linear function of the power of the 
internal power sources. 

In other words, the thermal acceleration of the Pioneer 
10 and 11 spacecraft due to heat anisotropically rejected 
off the vehicles can be modeled by a simple equation (see 
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FIG. 10: Modeling the exterior temperatures of a pair of 
RTGs using finite element software. Predicted values are 
within a few tenth of a degree from temperature readings 
obtained from flight telemetry. 



details in [III): 



but also incorporating other fictitious forces. This pro- 
gram could also be used to guide our investigation with 
JPL's Orbit Determination Program, similar to the ear- 
lier efforts in the study of the effect (see discussion in 
flii). 

This work is aimed to answer three questions: 

• Does such an approach yield values for the effi- 
ciency factors that are consistent with the values 
calculated using conventional methods? 

• Is it possible to obtain a good orbital solution incor- 
porating the thermal recoil force, without resorting 
to the use of other forces? 

• Is it possible to distinguish between solutions with 
a thermal recoil force vs. solutions that incorpo- 
rate other forces, such as a constant acceleration 
pointing towards the Sun? 

The results of this work to date, while not conclusive, 
are encouraging. Preliminary estimates of are consis- 
tent with the results obtained from finite-element ther- 
mal analyses, including the sophisticated on-going efforts 
at JPL and our earlier, simpler numerical integrations. It 
is our hope that the two approaches will complement each 
other, and that the accuracy necessary to determine the 
extent to which anisotropic thermal radiation from the 
spacecraft can be responsible for the anomalous acceler- 
ation (i.e., an acceleration accuracy of ^ 10^^° m/s^ will 
be achievable. 



«^thcrmal 



3mc ^ 



(3) 



where c is the speed of light, m is the spacecraft's mass. 
Pi is the i-th heat source, and is the associated dimen- 
sionless efficiency factor. (The additional factor of 2/3 
arises as a result of modeling the spacecraft's exterior 
surfaces as Lambertian emitters.) 

The values of Pi are known from flight telemetry or 
design documentation (Fig. [TT|) . The values of £,i can be 
determined approximately by analyzing the spacecraft's 
geometry and surface properties, as outlined in the pre- 
vious section. 

The simplicity of the relationship between heat and 
acceleration suggests that another approach may be pos- 
sible. Eq. ([3]) may be incorporated directly into the equa- 
tions of motion, allowing one to designate the efficiency 
factors £,i as parameters to be fitted. This, radically dif- 
ferent approach requires no knowledge of the geometry 
or thermal properties of the spacecraft. Question is, does 
it yield believable results? 

Recently, one of us (VTT) constructed an orbit de- 
termination program specifically to incorporate into the 
equations of motion Eq. ([3]), with the values of Pi sup- 
plied directly from flight telemetry. The software allows 
us to test a variety of hypotheses, attributing some, or 
all, of the anomalous acceleration to thermal radiation, 



V. CONCLUSIONS 

By 2007, the existence of the Pioneer anomaly is no 
longer in doubt. A steadily growing part of the commu- 
nity has concluded that the anomaly should be subject 
to further investigation and interpretation. Our contin- 
uing effort to process and analyze Pioneer radio-metric 
and telemetry data is part of a broader strategy (see dis- 
cussion in ^4, 5] ) . 

Based on the information obtained from flight teleme- 
try, we were able to develop a high accuracy thermal, 
electrical, and dynamical model of the Pioneer space- 
craft. This model will be used to further improve our 
understanding of the anomalous acceleration and espe- 
cially to study the contribution from the on-board ther- 
mal environment to the anomaly. 

It is clear that a thermal model for the Pioneer space- 
craft would have to account for all heat radiation pro- 
duced by the spacecraft. One can use telemetry informa- 
tion to accurately estimate the amount of heat produced 
by the spacecrafts' major components. The next step is 
to utilize this result along with information on the space- 
crafts' design to estimate the amount of heat radiated in 
various directions. 

This entails, on the one hand, an analysis of all avail- 
able radio-metric data, to characterize the anomalous ac- 
celeration beyond the periods that were examined in pre- 
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FIG. 11: Heat generated by RTGs (red, approximately 
straight line) and electrical equipment (green) in Pioneer 10 
over the lifetime of the spacecraft. 

vious studies. Telemetry, on the other hand, enables us 
to reconstruct a thermal, electrical, and propulsion sys- 
tem profile of the spacecraft. Soon, wc should be able to 
estimate effects due to on-board systematic acceleration 
sources, expressed as a function of telemetry readings. 
This provides a new and unique way to refine orbital 
predictions and may also lead to an unambiguous deter- 
mination of the origin of the Pioneer anomaly. 
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